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The major problem posed by the Fermi-Pasta-Ulam work was the relaxation (thcr- 
malization) of a many-body system. We review two approaches to this problem, the 
ergodic theory and Langevin (stochastic) equation ones, which have been applied 
for the description of statistical mechanics of stellar systems. Among the particular 
problems that we consider, are 1-dimensional long-range interacting systems and 
iterated maps and 3-dimensional N-body gravitating systems. 

1 The Fermi-Pasta-Ulam system: Historical introduction 

The numerical experiment by Fermi, Pasta and Ulam opened several new 
fields of research and appeared to be an invaluable source of inspiration for 
generations of physicists. 

This work is outstanding at least by the following reasons: 

1. It represents the first computer study of a nonlinear system; 

2. The results contradicted the belief held since Poincare. Fermi had con- 
sidered it 'a little discovery' (as quoted by Ulam), thus immediately evaluating 
their extraordinary importance; 

3. It was one of last Fermi's works, completed after his ultimate death in 
1954; 

4. Remained unpublished during a decade; 

5. Coincides in time with Kolmogorov's theorem (1954), though FPU and 
Kolmogorov-Arnold-Moser (KAM) theory 'met' only in 1966. 

6. Inspired the discovery of solitons and numerous other studies; 

7. Its results are not fully understood till now and the FPU model con- 
tinues its inspiring mission also today, after half a century. 

The circumstances around this work are now widely known, mainly due to 
the reminiscences of Fermi's collaborator of Los Alamos period and coauthor 
of this and other works, renown mathematician Stanislav Ulam. We mention 
some aspects only. 

The study was planned in the summer of 1952, during one of Fermi's 
visits to Los Alamos, the main calculations were carried on in the summer of 
1953, and final calculations in early 1955. It was documented as Los- Alamos 



gurzadyan: submitted to World Scientific on February 2, 2008 1 



internal report LA-1940 in May 1955 EJ. 

The calculations were performed on MANIAC, one of the first computers, 
together with the programming assistant John Pasta and with the help of 
Mary Tsingou. 

Though Ulam had reported the results in several meetings and had pub- 
lished their short account □, the report LA-1940 was fully published only in 
1965, in the second volume of collected papers of Fermi 111. After that, Ulam 
had studied other systems as well, i.e. described by homogeneous quadratic 
transformations, a class of so-called Binary Reaction Systems in three and four 
variables and of cubic transformations in three variables, in collaboration with 
P.Stein and M.Menzel (Tsingou) I. 

What was the aim of FPU work? 

According to Ulam □: 

"The ultimate aim was to discuss problems with more independent vari- 
ables in the hope of obtaining material which would suggest some general fea- 
tures of behavior of systems with an infinite number of degrees of freedom, with 
nonlinear interaction terms as they occur between the oscillators in quantum 
theory or between the degrees of freedom in an electromagnetic field or a meson 
field. The mathematical possibilities concerning what might be called quasi- 
states were discussed with Fermi in this connection" . 

FPU system represents a 1-dimensional system of N nonlinearly coupled 
oscillators given by a Hamiltonian 

N + l 1 1 

2=1 

at boundary conditions 

°j = Qj+N- 

The calculations were performed for three type of perturbations, cubic, 
quadratic and of broken-linear function. 

At m = 3 (a-FPU) and 4 (/3-FPU) and N = 16, 32, 64 with zero velocities 
of particles at t = 0, the system was run during time scales corresponding up 
to 80,000 cycles of the linear system. 

Instead of expected thermalization, i.e. a trend towards equipartition over 
modes, a regular behavior has been observed, namely, "mode 1 comes back to 
within one. percent of its intimal value so that the system seems to be almost 
periodic" U. 

Ulam recounts □: 

The results were entirely different qualitatively from what even Fermi, 
with his great knowledge of wave motions, had expected. 
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Indeed, it was contradicting the widely held views based on the Poincare 
theorem, that such perturbed systems have to be chaotic. 

By dramatic coincidence, at the same time thai belief was shaken by a 
theorem published without a proof by Kolmogorov in Moscow in 1954. Its 
proof was given almost a decade later, by Arnold and Moser in early 1960s, 
and Kolmogorov- Arnold-Moser (KAM) theory had to become an achievement 
of fundamental importance. 

Fermi was well familiar with Poincare theorem. Already in 1923 he had 
applied Poincare's method to prove that a system which has no other. integrals 
than the energy one, has to be ergodic on the energy hypersurface □. 

It might be, though Ulam mentioned about this indirectly p, that one of 
Fermi's motivations for FPU study was to check his theorem via an experi- 
ment. Does the fact that Fermi immediately accepted the paradoxical result of 
the experiment, indicate his anticipation of new fundamental insights? What 
would be his reaction, if informed about Kolmogorov's theorem? 

Unfortunately these questions remained unanswered. Fermi died, even 
before FPU was completed, without delivering his scheduled Gibbs Lecture 
for 1955 in American Mathematical Society. 

Kolmogorov's theorem needed around 15 years to reach the wider physical 
audience in west. 

Thus, FPU emerged both at the time of major breakthrough in the theory 
of perturbed nonlinear dynamical systems and at the moment of appearance 
of the first electronic computers. 

Fermi's intuition had clearly revealed both these aspects already at that 
time, as Ulam refers tl to Fermi's opinion on the importance of the 'under- 
standing of nonlinear systems ' for the future fundamental theories, and the 
'potentialities of the electronic computing machines ' and even mentions about 
Fermi's learning of the actual coding (programming) during one summer. 

The forthcoming decades indeed were marked by spectacular achieve- 
ments both in the studies of nonlinear systems and in the development of the 
computer technique. 

In 1958 Kolmogorov introduced the metric entropy, now known as 
Kolmogorov-Sinai (KS) entropy, as an invariant of the transformation, In 
1959 he introduced K-systems, now known as Kolmogorov systems^. 

Introducing the entropy as a metric invariant, Kolmogorov moved to prob- 
abilistic and geometric methods from the functional approach dominating 
since von Neumann and Hopf. 

His works together with those of Sinai and Anosov in 1960s are considered 



"Kolmogorov's original notation for K-systems implied 'quasi-regular systems'. 
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to lead to the modern ergodic theory of smooth systems! 

FPU work stimulated important studies, already before the publication 
in Fermi's volumes. We mention only two. Tuck and Menzel (Tsingou) in 
1961 observed a superperiod at more extensive studies of the same system. 
Their results remained unpublished until 1U72, though Ulam mentioned them 
when commenting Fermi's volume in 1965 u. StilLunpublished, FPU inspired 
the discovery of solitons by Zabusky and Kruskal Efl in 1965, "as solitary wave 
pulses which pass through each other and do not sum in a usual linear-wave 
superposition." 

In forthcoming years FPU continued and still is continuing to serae—as 
a reference point for numerous studies (for reviews and references seeB-Ej). 
Already the fact, that the decade delay of its publication, which might be 
fatal for any other scientific result, had not affected its influence, indicates its 
extraordinary power. 

In the present review, on the example of FPU, we will discuss the prob- 
lem of thcrmalization of nonlinear systems associated with an important astro- 
physical problem, the relaxation of N-body gravitating systems. The structure 
of the paper is as follows. 

First, we briefly discuss the association of FPU to the Poincare theorem 
and KAM theory. Then, we review the approach to the thermalization based 
on the Langevin equation, as of the simplest and most transparent version of 
the Brownian dynamics. Some results on the long-range interacting systems 
and iterated maps are presented thereafter and finally, we discuss the geomet- 
rical approach in stellar dynamics. Each topic is rather broad, and therefore 
we confine ourselves with brief accounts, without going into the details of 
applications. 

2 Relaxation in Statistical Mechanics 

Thus, the FPU puzzle concerned the thermalization, or relaxation towards 
equilibrium, one of the oldest problems in statistical physics. 

On the one hand, the relaxation verifies the equilibrium statistical me- 
chanics, since provides the basic equilibrium distribution functions — mi- 
crocanonical, canonical and grand canonical — as the result of a dynamical 
relaxation process. On the other hand, the problem of thermaliaa.tkm con- 
stitutes the very subject of non-equilibrium statistical mechanics Ej'EZI, since 
it is supposed to describe scenarios of relaxation and to give information on 
relaxation times scales. 

Historically, the first step in this field was made by Boltzmann. His ki- 
netic equation for weakly non-ideal dilute gases was the first physical example 
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which derived thermalization starting from the microscopies dynamicsQ'0. 
In Boltzmann's setup one considers a closed macroscopic system, and makes 
certain sensible kinetic assumptions on its dynamics which result in the ki- 
netic equation. This is sufficient to derive the H-theorem, which states a 
monotonous increase of entropy, and as the main result one gets relaxation of 
the macroscopic system towards equilibrium (thermalization). 

Unfortunately, this approach is not always easy to motivate and imple- 
ment. In particular, the kinetic assumptions become mathematically rigor- 
ous only in rather artificial limits. For instance, it is known that a self- 
consistent description of thermalization for essentially non-ideal gases within 
Boltzmann's approach meet serious problems 113. 

The example of Fermi-Pasta-Ulam comes to verify that not all realistic 
systems need to satisfy conditions of Boltzmann's kinetic equation. 

Another approach to the thermalization problem was started by Einstein 
and Langevin EaOEJ. This is a drastic change in philosophy as compared 
with the Boltzmann's approach, since one does not look for thermalization 
of the whole macroscopic system, but rather investigates its subsystem. This 
approach originated from the problem of Brownian motion, where one follows 
to the motion of one heavy particle immersed in a thermal bath. 

In contrast to the Boltzmann's one, in this approach the thermalization 
problem appears to be more transparent and more widely applicable, since 
after certain reasonable and rather universal assumptions on the structure of 
the bath, one derives, theijrnalization of the Brownian particle via an exact 
dynamical analysis ll3'EZHl3ll3. 

As a consequence, one gets relaxation towards Gibbs distribution. One of 
the essential advantages of this approach is that the extension to quantum- 
mechanical systems appears to be straightforward lii 

The third approach is provided by ergodic theory which classifies dynam- 
ical systems by so called statistical properties, including those of tending to 
the equilibrium state. 

In particular, the property of mixing is defined via decay of time corre- 
lation functions, as necessary condition for relaxation, while the ergodicity 
alone is not sufficient for that. 

Though this approach is known by difficulties of checking the sufficiency 
conditions for realistic physical problems, there are methods, among them are 
the geometrical ones, which enable to overpass those difficulties. 

The geometric methods of ergodic theoaz. have been advanced by Kol- 
mogorov, Sinai, Anosov and others_(see EEl), while in physical problems 
they were firstly applied by Krylov E3 based on the early results of Hedlund 
and Hopf. 
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In stellar dynamics, i.e. for the problem of N-body gravitating systems 
interacting with Newtonian potential, the Langevin equation approach has 
been used long time ago by Chandrasekhar and von Neumann The 
geometrical methods ,p£-ergodic theory firstly have been applied in stellar 
dynamics in O[^0c3E3e3 (and their other papers quoted therein). 

3 Poincare theorem, KAM and FPU 

First, let us recall that N-dimensional system is considered as integrable if its 
first integrals 7i,...,//v i n involution are known, i.e. their Poisson brackets 
are zero. Then, as follows from the Liouville theorem if the set of levels 

M I = {I j (x) = I$,j = l,...,N} 

is compact and connected, then it is diffeomorphic to N-dimensional torus 

T N = {(6i,...,6 N ),modd27r}, 

and the Hamiltonian system performs a conditional-periodic motion on Mj. 
Poincare theorem El states that for a system with perturbed Hamiltonian 

H(I,<p,e) = H a {I) + eH 1 (I,<p,e), (2) 

where I, tp are action-angle coordinates, at small e > no other integral ex- 
ists besides the one of energy H — const, if Hq fulfills the nondegeneracy 
condition, 

det\du/dl\ + 0, (3) 

i.e. the functional independence of the frequencies lo — dH^/dl of the torus 
over which the conditional-periodic winding is performed. 

Though this theorem does not specify the behavior of the trajectories of 
the system on the energy hypersurface, up to 1950s it was widely believed 
that such perturbed systems have to be chaotic. 

Kolmogorov's theorems of 1954, now the main theorem of KAM theory, 
came to contradict to that belief, namely, formulating the conditions when 
the perturbed Hamiltonian systems can remain stable. 

It states: 

If the system (2) satisfies the nondegeneracy condition (3) and H± is an 
analytic function, then at enough small e > most of non-resonant tori, i.e. 
tori with rationally independent frequencies satisfying the condition 

5> Wfe ^0, (4) 
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do not disappear and the measure of the complement of their union set 
H(M) -> at e -> 0. 

The theorem says nothing about the limiting value of the perturbation 
e. Later numerical studies showed that, for certain systems that theorem 
can remain valid even at not small perturbations. Also, the behavior in two- 
dimensional Hamiltonian systems when the torus separate the regions of three- 
dimensional system and those of higher dimensions have to be different, when 
the phase trajectories of the system can wander to other regions due to Arnold 
diffusion. _ 

Izrailev and Chirikov in 1966 113, were the first to match KAM and FPU, 
via the formulated resonance overlap (Chirikov) criterion. Though even ear- 
lier, in 1963, Ford and Waters quoted Kolmogorov's address at the Interna- 
tional Congress of Mathematicians of 1957: 

''Kolmogorov claims to have shown that nonlinear systems, more general 
than the systems considered here, are not ergodic", the importance of that 
theorem was not been recognized and the works of Arnold and Moser were 
unknown. 

Now, after numerous studies it is clear that FPU does not satisfy the con- 
ditions of the KAM theorem, though its behavior is close to that of integrable 
systems and to KAM conditions. KAM-FPU association is still an open issue. 

FPU does not satisfy Poincare theorem conditions either, and cannot 
be directly traced via Poincare recurrence time, as was attempted by some 
authors (see refs in 13) . 

The rigorous understanding of the thermalization remains the challenge 
of FPU. 



4 The Langevin equation 

Let us first consider the derivation of Langevin's equation within the Brownian 
motion (system-bath) approach. Instead of describing this derivation in its 
full generality, we will consider the simplest situation as to provide all essential 
mathematical details, and to make clear all necessary physical assumptions. 

The Langevin equation is derived from the exact Hamiltonian description 
of a subsystem (Brownian particle) and a thermal bath, by tracing out the 
degrees of freedom of the bath. The influence of the particle on the bath is 
assumed to be sufficiently small. Thus, only the linear modes of the bath are 
excited, and the interaction of the particle with the bath is assumed to be 
linear. To be as pedagogic as possible, we first take a definite model for the 
bath, namely a collection of harmonic oscillators. 
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For the total Hamiltonian we thus assume 

P 2 



Htot = H + Hb + Hi = jj- h V(x) + V 

2m ' 



2m,- 2 * 



E 



(5) 



where if, Tie are the Hamiltonians of the particle and the bath, and Hj is 
the interaction Hamiltonian. p, pi, x, Xi are the momenta and coordinate 
of the Brownian particle and the linear modes of the bath. Ci are coupling 
constants with the bath; their concrete expressions will be given later on. V{x) 
is the confining potential of the particle, and m and rrii are the corresponding 
masses. 

The Brownian particle was taken to be one-dimensional. This is by no 
means a restriction of generality, because the generalization towards three- 
dimensional particle, as well as N interacting Brownian particles, is straight- 
forward. 

Notice that the a^-terms form a complete square, since it includes a self- 
interaction term proportional to x 2 . This guarantees that the total Hamilto- 
nian H t ot will be positive definite. 

At the initial time t = the bath and the particle are decoupled. The 
distribution function of the particle is not specified. 

The bath, on the other hand, is assumed to be in equilibrium at tempera- 
ture T — 1/(3 (we set Boltzmann's constant ks = 1 in the body of this work), 
and described by its Gibbs distribution. 

This means 



Psb(0) = p(0) PB (P) = ps(P) 



cxp[- PH B (0)] 
Z B 



(G) 



where p(0), ps(0) are the initial distribution functions of the particle and the 
bath, and 



Z B = / II dxi exp[-f3H B ] 



is the partition sum of bath. 

The Hamilton equations of motion for the bath modes read 

. _ 1 

Xi Pi , 



m, 



-XiUliWi + CiX 



(7) 
(8) 



gurzadyan: submitted to World Scientific on February 2, 2008 8 



These equations are solved readily: 

Xi(t) = Xi(Q) cosujit H — sinwitH — / ds smu>i(t — s)x(s) (9) 

»'.-', »',-', ./[• 

Pi(t) — — rriiU)iXi(0) siuuJit + pi(0) cosujit + Ci / ds cos u>i(t — s)x(s) (10) 

Jo 

Thus, the Heisenberg equations of motion for the Brownian particle read 

x=-p, (11) 
m 

j> = -V'{x) + Y i c i x i -xY i —± s (12) 

i i L 

Combined with Eq. (^) the last equation becomes 

p=-V'(x)-x(0)j(t)- [ dsj(t- s)x(s) + r)(t), (13) 



n 



where 



rj(t) — Ci[xi(0) cosujit + ^ sinujit] 
t— 1 rriiUJi 



E \/ ^T^ ^ + ^(0)^^*], (14) 



7W = E^ cos (^)' ( 15 ) 

are the noise related to the unperturbed bath, and the friction kernel, respec- 
tively. Notice that in this exact derivation the back-reaction of the bath on 
the particle has been taken into account. 

It is described by the integrals in eqs. (||), (|l(i|), and brings the damping 
terms x(t)j{0) - x(0)-/{t) - f* dsj(t - s)x(s) in cq. ([[J) to yield (|l|). 

4-1 Drude-Ullersma spectrum 

For some, but not all, applications it is benefitable to consider a fully explicit 
case for the bath. The bath is assumed to have uniformly spaced modes 

u>i=iA i = 1,2,3,- •■ (16) 
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and for the couplings we choose the Drude-Ullersma spectrum 



C ^V^^^TT^ (17) 

Here T is the characteristic Debye cutoff frequency of the bath, and 7 stands 
for the coupling constant; it has dimension kg/ s. 

The thermodynamic limit for the bath is taken by sending A — > 0. Notice 
that this creates an infinite "Poincare" timescale 1/A, implying that in the re- 
maining approach the limit of "large times" always means the quasi-stationary 
non-equilibrium state where time is still much less than 1/A. 

In the limit A — > each coupling Ci ~ \fK is very weak. The fact that 
the bath has many modes nevertheless induces its non-trivial influence. 

At finite but small A the system has an initial relaxational behavior, 
which at times of order 1/A is changed in a chaotic behavior. 

It is customary to define the spectral density 



UliLJi 



It has the Ohmic behavior J w juj for lu <C T, and 7 is called the interaction 
strength or damping constant. 

As J(io) is cut off at the "Debye" frequency T, it is called a quasi-Ohmic 
spectrum. 

It can then be shown that the friction kernel (|l5|) becomes 

2-7 f°° r 2 

7(0 = — / dco cos^^7re- r ' t ' (19) 

It is non-local in time, but on timescales much larger than 1/r it may be 
replaced by j5 + (t). 

The resulting Langevin equation reads 

p + — f dt'e-^-^plt') + V'ix) = - 7 re- rt x(0) + r)(t), (20) 
m Jo 

It follows from @ that the noise is stationary Gaussian, and has the following 
properties: 

K { t-t>) = {nm t>)) = H f^^| 7f §=27T§e-l-'l (21) 

n j 1 + (uj/ry 2 

where, again, the average is taken over the initial state of the bath (over 
realizations of the noise). 
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The connection between properties of the noise^pJ-U^fxiction kernel is 
the consequence of fluctuation-dissipation theorem Oli^O'E^lli 

The classical white noise situation Kit) — *■ 2~fT5(t) is recovered by taking 
the high-r limit (Ohmic spectrum of the bath frequencies). 

In this Ohmic regime one gets finally for t > 0: 

P + ^P(t) + V'{x) = rj(t), (v(tMt')) = *yT8(t ~ t'). (22) 

Finally we wish to, rnontinn that there are alternative ways to derive the 
Langevin equation II3Ea, since many of its properties are rigidly determined 
by general statements like the fluctuation-dissipation theorem Ea. 

Nevertheless, we choose to focus on concrete models, because they show 
in detail how the Langevin equation arises from first principles, and thus are 
better suited for pedagogical purposes. 

4-2 Fokker-Planck equation 

Once we are interested by one-time quantities, the state of the Brownian 
particle is described by the distribution function p(p, x, t): 

p(p,x,t)= J dp dxop(p ,Xo,0){8(p{t)-p)5(x{t)-x)), (23) 

where the average is taken with respect to the noise, p(p, x, t) and pipo, xq, 0) 
are final and initial distribution functions, while p{t), x(t) are the solutions of 
( p2| ) for the corresponding initial conditions, and for a particular realization 
of the Gaussian noise. 

The Gaussian noise is distributed according to the functional 

V[rf\ ~ exp-^y dtdsr)(t)K- x (t - s)rj(s). (24) 

It will be clear that the distribution function W(yi, yi,t), and the propagator 
(6(p(t) — yi)S(x(t) — 2/2)) satisfy the same differential equation although with 
different initial conditions. 

Therefore, we will use them interchangeably, making additional specifica- 
tions only when it will be necessary. 

Differentiating p[yx,y2,t) we get 

dpjyi,V2,t) \- djv k p) d 

— dl — = - L -q^~ ~ dy-^ 5 ^ ~ - i*M*)>. (25) 



where 



Vl = -l p (t)-V'(x), v 2 = ^-. (26) 
m m 



gurzadyan: submitted to World Scientific on February 2, 2008 11 



are the damped Newtonian acceleration and the velocity, respectively. 
The following two relations can be verified 



v(t)V[ V }=- / dt'K(t'-t) 



(27) 



Sp(t)/Sr,(f) 



9(t-t') jexpj^ duA(x(u))j 



(28) 



A(x) = 



-7/m -V"(x) 
1/m 



(29) 



where {■■■}+ means the chronological ordering (time-ordering). 

The first equality holds since V[rj] is Gaussian; the second follows directly 
from Eq. (E2) because of the relations 



d 6p(t) 7 Sp(t) Sx(t) , 

dtd V (t')~ m5r)(t') 1 %(f) 
d 5x(t) 1 5p(i) 



di 5r}{t') m Sr)(f) 
Using Eqs. (|5|, |7|, ||, ||) we obtain 

dp(yi,V2,t) 



(30) 
(31) 



<9t 



= + ^ - viMx(t) i&)<M{*}, *)>} , 



where $ is the 2x2 matrix 



$({x},i) = y di'iT(t') |exp^ duA(x(u)) 



(32) 
(33) 



and is the corresponding matrix element K(t') = 2~fT5(t'), implies 

*ki=-yT8 kl . (34) 

is thatp-wfi-iObtain a diffusion-type equation (Kramers-Klein 
itself HW: 



The final result is t 
equation) for p 



^£^lEl3. = -L2E + —([l p + v'(x)]p)+'yT— (35) 
dt m dx dp V m ) dp 2 
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4-3 Thermalization and H-theorem 



It is now possible to provide a general argument that Eq. ( |35| ) ensures ther- 
malization. First one notices that its stationary solution, obtained from ( |35| ) 
by setting dtp = 0, is given by the standard Gibbs distribution: 



Pst 



|exp 



1 2m 



where Z is the statistical sum. Second, one notices that the H-function 



H(t) — J dpdxp(p,x,t) In 



Pip, x, t) 

p st ip,X,t)' 



(36) 

Emm 



(37) 



where pit) is the actual time-dependent solution, monotonously decreases to- 
wards zero, which is its equilibrium value p = p s t- 

Indeed, differentiating Hit), using ( f35| ) and requiring that both p and p st 
are zero at infinity, one gets: 



Hit) — ~jT / dpdx p{p 7 x,t) 



d p(p,x,t) 



n 2 



dp p st (p,x,t) 



< 0. 



(38) 



Thus, all solutions of 
bution (H 



J) relax with time towards the stationary Gibbs distri- 



Eq. (38) provides an information on the global relaxation time; 
it is proportional to the damping constant 7 and temperature T. Notice that 
the 

H-function provides an. iiifarmation-theoretic distance between the distri- 
bution functions p s t and p 

Thus, the open system approach provides a rigorous justification of ther- 
malization. Let us outline the main assumptions of the present setup. 

1) The Brownian particle was assumed to be in contact with a macroscopic 
thermal bath at equilibrium. 

2) The influence of the particle to the bath was assumed to be weak 
enough, so that the non-linear modes of the bath are not excited. 

3) The choice of initial conditions (|^) was special, and moreover, it has 
to be special. One can motivate that for any macroscopic, system there are 
initial states which do display anti-relaxational behavior llJj—Por sufficiently 
simple cases one can show that those states are exceptional E3. 

4) The thermalization was obtained via a rigorous approach; no a priori 
assumptions on ergodicity and mixing were used. As the main result Eq. ( |38| ) 
provides a general expression for the relaxation time. 
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5 Long-Range Interacting Systems 

The problem of understanding of the observed thermalization in FPU system 
thus includes the explanation of the existence of large relaxation times tr 
and their dependence on the number of particles N and the energy density 
e = E/N (see 113 and references therein). 

We will mention some of the remarkabLa-results obtained in the recent 
period. The following dependence was foundL3 



TROL^f, (39) 

for finite ./V and fixed wave number interval for the /3-FPU model, i.e. the 
divergence of the equipartition time scale with increase of N. On the_pther 
hand, in the thermodynamical limit the dependence on N disappears E3 

tr cx e- 3 . (40) 

The conclusion of these studies is that, "the process of relaxation to equiparti- 
tion in the FPU model is not regulated by the microscopic chaotic instability 
but by the typical time in which an orbit diffuses iu phase space, which is 
determined by the interaction among the phonons" t3. 

While moving from FPU to 1 and 2-dimensional models with long-range 
interaction, to mimic the gravitational systems, the problem becomes more 
complex since acquires new thermodynamics, negative specific heat, etc. 

Indeed, different results have been obtained at the analysis of the chaotic 
properties on the following 1-dimensional 



TV 



2 



h = T,^ + JnE n ( 1 - cos (^ - ( 41 ) 

i—1 — 1 



and 2-dimensional systems 



^ 2m y ' 

i—1 

+ N $ ~ C0S ( Xi ~ x i) ~ cos (^ _ ^) ~ cos ( x i ~ Xj)cos(yi - yj)], 

and, particularly of their equipartition time scales. As a result, "few scaling 
laws have been found to be universal for all these models. In the gaseous 
phase the maximal Lyapunov exponent vanishes as 

iV~ 1/3 (43) 
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in analogy found inE13 for some specific gravitational systei 
In a similar model, the Konishi-Kaneko iterated map t 



N 

fc^sin27r(a;?-0, (44) 
i=i 



x ™ +1 =x?+p? +1 ;(modl). (45) 
which can be represented as a system with Hamiltonian 

H = E N lY + ^ E Ncoapirfa Xj ))] (46) 
i=i ijtj 

the appearance of cluster- type configurations, of power law conflations El, and 
Feigenbaum period-doubling bifurcations have been observed cB. Particularly, 
the bifurcation scale 6 = 8.72... can be estimated at the following conditions 
for the bifurcation points 

iv 

y £\xp 1 -x?\<e,(2 1 =2), 

3 



N N 

£ \x] +2 - x]\ < e^T" ~ < e ' ( 22 = 4 )' 

3 3 

for each k n , n = 1, 2, respectively, where e is the accuracy of the obtained 
values of k n . 

The use of iterated maps can enable one to avoid the principal difficulties 
associated with iV-body systems - the non-compactness of the phase space 
and singularity of Newtonian interaction - which are not always crucial at the 
study of a particular property of a given stellar system. Also, as we saw for 
Konishi-Kaneko map, the iterated maps can be informative in revealing the 
mechanisms of developing of chaos. 



6 Gravitating Systems 

While moving to the gravitating systems first one must recall that, the Hcnon- 
Heiles lid system being one the first demonstrations of chaos at numerical 
experiments, was motivated by the study of the motion in the axisymmetric 
Galactic potential. 



gurzadyan: submitted to World Scientific on February 2, 2008 15 



Since then, developments in nonlinear dynamics had various impacts pn 
|tl|e problem of N-body gravitating systems (for reviews see the volumes ej, 

The chaotic dynamics had offered principally new insight on-the properties 
and the pplution of the Solar system. In the works of Laskar liJ. TremaineEJ, 
Wisdomc3 and their collaborators it was shown that the inner planets can 
possess chaotic properties influenced by the perturbations of the large planets. 
Particularly, Mercury and Venus undergo chaotic variations of the eccentricity 
of their orbit, which can lead to overlapping of orbits with subsequent drastic 
change of the orbit of Mercury, up to its escape from the Solar system. 

By means of the frequency map analysis technique developed by Laskar, 
the stabilizing role of the Moon in the chaotic variations of the obliquity of 
the Earth has been shown with direct consequences for the climate and hence 
for the evolution of life on the Earth. Chaotic variations of the obliquity were 
shown to be important for Mars. ._. 

The idea of the frequency map technique Cj is in the search of quasiperi- 
odic approximations of the solutions of the perturbation problem by means 
of finite number of terms 

N 

Xj(t) = x" + a(m k ; v) exp[i < m k ,v > t] (47) 
fe=l 

where Vj are the frequency vectors. Note, that the system is close to an 
integrable one, and possess KAM (invariant) tori in the phase space. 

So, proper theoretical treatment on the perturbed systems together with 
advanced computations had led to such basic results. 

While profound role of chaotic effects in the evolution of Solar system is 
established, their role in stellar dynamics is often remained underestimated^] 
even though the latter ones are not close to integrability and KAM conditions, 
and one might expect that chaos has to be of more importance for them than 
for planetary systems. 

The problem of relaxation of stellar systems can illustrate the situation. 

The two-body relaxation formula derived many decades ago via the linear 
sum of contributions of several two-body interactions, is still the basic refer- 
ence for stellar system, even though it predicts time scales exceeding the age 
of such well-relaxed systems as the elliptical galaxies. 

A priori it is not excluded that the two-body relaxation time scale can be 
the solution of the nonlinear problem even when the perturbations of other 

'See e.g. 0, where the authors needed half a page to transfer from linear to nonlinear 
effects in describing the relaxation of stellar systems. 
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particles are included, however this does not follow from the 'linear' deriva- 
tion borrowed from plasma physics, involving the Debye screening and the 
Coulomb logarithm, i.e. ignoring the long range interaction of gravitating 
particles. 

While even for 1-dimensional systems the relaxation is essentially different 
for short and long range interactions. 

Geometrical methods appeajLaaa useful tool for the study of nonlinear ef- 
fects of gravitating systems. In EjJHj the K-mixing and exponential tending to 
an equilibrium state (quasi-equilibrium) was shown analytically for perturbed 
spherical N-body systems. 

That was reached via the proof that while the two-dimensional curvature 
of the configuration space for N-body systems is sign-indefinite, for spherical 
systems in large N limit it is strongly negative. 

The formula derived for the relaxation time scale for the parameters of 
real stellar systems differs from the two-body relaxation formula, and by now 
is supported by (a) alternative theoretical derivation as mentioned above 
lb) numerical simulations and (c) observational data of star clusters (see I 
Ej and references therein). 

The importance.^ chaotic effects in relaxation of stellar systems was 
shown by PfennigeiO using the Lyapunov exponent technique. That study is 
especially remarkable given the. limited applicability of Lyapunov exponents 
for many-dimensional systemsEa. I— . 

Back in 1964 at a remarkable computer study of his time Miller c3 had 
observed exponential growth of errors of the coordinates and velocities, A = 
J2(^x 2 + Sv 2 ) for systems with up to 32 gravitating particles. 

The result of Miller remained not well understood for quite a long time. 
Recently, in some papers (e.g.E3 and refs therein) it is identified with the ex- 
ponential instability (K-mixing) shown in even though their exponential 
time-scales are quite different. Since such opinion is moving from one paper 
to another, we will consider this issue in more details. 

As follows from the Maupertuis principle, only trajectories in the region 
of the pOpjELfiguration space M = [E — U(r) > 0] defined by the Riemannian 
metric! 



50 



N 3 

ds 2 = {E- f/(r M ..., r N , 3 ) £ £ ( 4 ») 

a— 1 z—1 

are corresponding to the N-body Hamiltonian system 

N 3 2 
a—1 i—1 
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s m 
for 



where 

„ M \- Gm a m b ) 

a<b ' 1 

The Maupertuis parameterization and hence Riemannian metric guarantee 
the conservation of total energy of the system 

H = E = const, 

and therefore enables to study the hyperbolicity and mixing of the system 
with exponentially deviating trajectories. 

The works discussing the exponential instability can be divided into two 
groups: p 

1. —Based on Maupertuis paramptpr.ization: by Krylov, AnosovEl and 
ArnoldEa, and in stellar dynamics, Ej'ESEa; 

2. Without using Maupertuis parameterization: numerical experiments 
by MillerEil, repeated by Lecar and Standish in 1960s and again in 
many papers by Kandrup, Goodman, Heggie, Hut and others (see 
references) . 

In the first group of works the exponent is linked with mixing and re- 
laxation via theorems of ergodic theory. For the second group, no such link 
exists. The reason is as follows. 

Without Maupertuis parameterization, two points on two phase trajecto- 
ries can deviate while the trajectories themselves may not deviate at all. Such 
apparent exponential deviation will be observed for systems of any number 
of particles, even of very few, and of amz_configurations, say of disk type. 
In contrary, the exponential deviation ofOtJ exists only for spherical systems 
and only at large N-limit, and not for disk or few-particle gravitating systems. 

Therefore, a statement that, the second group papers also refer to chaos 
and mixing, is the same as to assign chaos to the two-body system, where the 
same apparent deviation can be observed. 

This is the reason that, only the dynamical (crossing) time scale is ob- 
served at such numerical experiments. 

Numerical simulations taking into account the Maupertuis parameter- 
ization are in agreement with the relaxation time scale derived in Ea€a. 

At the same time, we are not aware of any numerical experiments which 
confirm the two-body relaxation time scale^. 



c Donald Lynden-Bell to whom we consulted (V.G., July,2001), was not aware of such ex- 
periments either. 
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The fact that Miller's experiment does not imply mixing, by no means 
diminishes the pioneering character of his work. However, now the results of 
theory of dynamical systems provide its interpretation. 

The geometrical approach enabled to arrive also to other general conclu- 



elliptical and spiral galaxies should have different origins. 

The spiral galaxies particularly. are more regular than ellipticals, i.e. do 
not possess the property mixing E3. Due to the existence of fundamental 
frequencies they can be described by Laskar's frequency map technique E3 

At Geneva workshop in 19113 one of us (V.G.) presented a list of "10 key 
problems of stellar dynamics"EJ. Now, about 10 years after, one can note 
that though a considerable work has been done, especially a numerical one 
due to advances of computer technique E3, still there is a long way to go to 
fully evaluate the role of nonlinear effects for stellar systems. 

7 Epilogue 

The Fermi-Pasta- Ulam problem was able to address profound issues of nonlin- 
ear systems, though studied mainly from physical concepts, rather than from 
rigorous mathematical ones. It still lacks answers to basic questions. As men- 
tioned Pesin E3: " It is a very interesting question the relation between KAM 
and FPU and whether FPU can serve as an " experiment" for KAM" . One 
can expect that the activity in FPU problem will shift towards investigation 
of its rigorous aspects. 

The relaxation problem remains a central issue for many systems, includ- 
ing open and long-range interacting ones. In astrophysical context it concerns 
N-body gravitating systems, galaxies and star clusters. We briefly discussed 
two powerful approaches which have been applied in stellar dynamics, the 
Brownian particle approach and the geometrical methods of ergodic theory. 

In stellar dynamics also one has to expect more rigorous studies of non- 
linear effects, including in the performance and interpretation of numerical 
experiments, a lesson given to us by Fermi, Pasta and Ulam 50 years ago. 
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